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Automatic  Detection  of  Hail  by  Radar 


I  INTRODUCTION 

There  has  long  been  a  need  lor  an  operationally  effective  means  of  distinguishing 
storms  that  produce  hail.  Over  the  vears.  radar  meteorologists  have  been  able  to 
identify  hail -producing  cells  with  some  success  bv  examining  radar  returns  for  one 
or  more  storm  characteristics  such  as  echo  tops,  maximum  reflectivities,  heights 
of  various  reflectivities  and  tilt.  Recently,  Lemon  performed  a  study  of  the  exist¬ 
ing  hail  identification  techniques;  he  derived  a  set  of  verv  successful  severe  hail 
storm  identification  criteria  lor  the  WSR-57  radar  based  on  the  three-dimensional 
reflectivity  structure  of  a  model  severe  hailstorm. 

The  structural  nature  of  the  Lemon  hail  criteria  suggested  the  development  of 
a  hail  analysis  algorithm  for  automated  identification  of  hailstorms.  The  types  of 

data  required  to  test  the  conditions  in  the  criteria  are  readily  available  in  our 

2 

existing  Modular  Radar  Analysis  System  (MRAS).  In  the  course  of  development, 
we  added  six  additional  hail  indicators  to  supplement  the  Lemon  criteria.  Four  of 
these  proved  useful. 

(Received  for  publication  29  September  1982) 

1.  Lemon,  L.  R.  (1978)  On  the  Use  of  Storm  Structure  for  Hail  Identification, 

Preprints,  18th  Conference  on  Radar  Meteorol. ,  Boston,  Am.  Meteorol.  Soc. , 
203-206. 

2.  Forsyth,  D.  E. ,  Bjerkaas,  C.  L. ,  and  Petrocchi,  P.J.  (1981)  Modular  Radar 

Analysis  Software  System  (MRASS),  Preprints,  20th  Conference  on  Radar 
Meteorol.,  Boston,  Am.  Meteorol.  Soc.,  696-6^7! 
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2.  HAIL  INDICATORS 


Our  nine  hail  indicators  were  formulated  to  determine  the  existence  of  radar 
reflectivity  structural  features  generally  associated  with  hail  producing  cells. 
These  hail  features  are  represented  schematically  in  a  model  of  an  ideal  hail  cell 
(Figures  1  and  2). 


Figure  1.  Vertical  Cross  Section  of  Radar 
Reflectivities  in  a  Model  Hailstorm 


Descriptions  of  the  nine  hail  indicators  and  the  rationale  for  their  use  are 
presented  as  follows: 

1-  Mid-level  (5  to  12  km  reflectivity  is  at  least  50  dBZ.  Indicators  1  through  3 
are  almost  direct  adaptations  of  the  Lemon  criteria.  Indicator  1  tests  for  a  given 
reflectivity  value  in  the  mid-level  regions.  Lemon  found  from  studies  by  Donaldson'^ 


3.  Donaldson,  R.  J.  (1961)  Radar  reflectivity  profiles  in  thunderstorms, 
J.  Meteorol.  18:292-305. 

-  W« 
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Figure  2.  Hadar  Slice  Schematic  of  the  Model  Hailstorm 
Shown  in  Figure  1 

and  Dennis  et  al  *  that  the  height  of  a  given  reflectivity  is  more  significant  than  the 
oeak  reflectivity  values  for  indicating  the  occurrence  of  hail.  We  found  this  to  he 
so  in  our  preliminary  testing,  but  we  increased  the  threshold  value  suggested  bv 
Lemon  from  45  dBZ.  to  50  dBZ.  to  decrease  the  number  of  false  alarms. 

2.  Mid-level  30-dBZ  contour  extends  at  least  4  km  beyond  the  40-dB/.  lowest 
level  contour.  Indicator  2  tests  for  the  overhang  of  the  mid-level  echo  over  the 
edge  of  the  low-level  echo.  This  test  checks  for  the  existence  of  a  weak  echo  region 
(WER)  indicative  of  strong  updraTts  usually  associated  with  severe  hailstorms.  * 

The  WKR  is  shown  in  Figure  1  as  a  region  of  reflectivity  loss  than  ,10  dBZ  located 
under  a  mid-level  reflectivity  contour  that  is  at  least  30  dll/..  As  a  result  of  our 
preliminary  tests,  we  altered  the  extent  of  the  overhang  in  the  Lemon  criteria  from 
6  km  to  4  km. 

3.  The  echo  top  lies  over  a  mid-level  overhang.  Indicator  3  tests  for  the  posi¬ 
tion  of  the  echo  top  relative  to  the  low  level  echo  |x>sition.  While  indicator  2  tests 
lor  the  existence  of  a  persistent  updraft  at  low  levels,  indicator  3  tests  for  its 
relative  strength  at  upper  levels.  Lemon*  concludes  that  if  we  have  a  strong 
updraft  in  a  strongly  sheared  environment,  we  should  find  the  top  above  the  weak 


4.  Dennir,  A.S.,  Smith.  1*.  L.  ,  .Jr.,  Ilovd.  ILL.  and  Musil.  U..I.  (1071)  Hadar 
Observations  of  Hailstorms  in  Western  Nebraska.  South  Dakota  School  of 
Mines  and  Technology,  Final  Heport,  NfS  FT!  rant"  t;  A -15  lft,  4  2  op. 


echo  region  or  above  the  overhang  itself.  These  conditions  are  shown  in  our  hail¬ 
storm  model. 

4.  The  maximum  reflectivity  at  any  level  is  at  least  55  dBZ.  Indicator  4 
searches  for  a  peak  reflectivity  value  above  a  given  threshold.  This  type  of 
indicator  has  been  used  in  the  past  with  varying  degrees  of  success. 

5.  The  area  of  the  30-dBZ  echo  at  any  mid-level  is  at  least  10  percent  greater 
than  the  area  of  the  lowest  level  30 -dBZ  echo.  It  was  the  intent  of  indicator  5  to 
supplement  indicator  2.  From  our  hail  model,  it  would  appear  that  the  area  of  a 
mid-level  slice  would  always  be  larger  than  the  area  of  a  low-level  slice. 

6.  The  direction  of  tilt  lies  to  the  right  of  and/or  behind  the  direction  toward 
which  the  cell  is  moving.  The  direction  of  tilt  is  defined  as  the  horizontal  angle 
between  the  direction  of  the  cell  motion  and  a  line  directed  through  the  centroid  of 
a  mid-level  slice.  A  centroid  is  'he  calculated  center  of  mass  for  a  given  slice  as 
shown  bv  lower  case  letters  a,  b,  d,  and  e  in  the  hail  model.  Indicator  6  supple¬ 
ments  the  test  bv  indicator  3  and  confirms  the  existence  of  a  WER  on  the  right  or 
rear  flank  of  a  cell,  as  shown  hv  the  model  hailstorm. 

7.  The  30-dBZ  top  is  at  least  8  km.  Echo  tops  have  long  been  used  as  a  means 
if  detecting  hail.  One  problem  is  the  selection  of  proper  height  thresholds  which 

varv  in  different  geographic  regions  and  with  the  height  of  the  tropopause.  The 
combination  of  reflectivity  threshold  and  height  (30  dBZ,  8  km)  was  selected  for  the 
Oklahoma  storms  on  the  basis  of  preliminary  tests. 

8.  The  direction  of  tilt  is  towards  the  south.  This  test  supplements  the  test 
hv  indicator  6.  If  we  assume  that  most  cells  move  in  an  easterly  direction,  this 
indicator  tests  for  a  cell  tilt  to  the  right  of  the  cell  direction  as  does  indicator  6. 
Indicator  8  covers  cases  where  the  direction  of  movement  for  a  given  cell  is 
unknown. 

9.  The  movement  of  a  cell  is  to  the  right  or  left  of  the  mean  motion  of  all  cells. 
Indicator  9  was  introduced  as  an  experimental  indicator  to  determine  if  hailstorms 
have  a  preferred  direction  of  motion  relative  to  the  direction  of  motion  of 

>ther  storm  cells. 


3.  ALGORITHM  TESTS 

The  hail  algorithm  is  one  of  four  Special  Analysis  Modules  (SAMS)  operating 

2 

under  the  Modular  Radar  Analysis  System  (MRAS).  “  I" he  MRAS  also  preprocesses 
the  radar  data  acquired  during  an  antenna  volumetric  scan  sequence  and  saves  these 
data  for  use  by  the  SAMs.  The  volumetric  scan  consists  of  a  predetermined  number 
>f  azimuthal  scans  at  selected  elevation  steps.  A  typical  volumetric  scan  sequence 
that  we  used  for  our  data  acquisition  was  composed  of  five  360°  azimuthal  scans  at 
elevations  of  0.  2,  I).  5,  1.  5,  3.  0,  and  5.  0  degrees. 
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Representations  of  tae  data  acquired  by  the  '-olumetric  scan  sequence  are 
shown  in  Figures  1  and  2.  The  straight  lines  represent  elevation  radials  that 
generate  slices  through  a  given  cell.  For  each  slice,  the  MRAS  provides  data 
defining  tne  slice  elevation  angle,  cross  sectional  area  enclosed  bv  tne  30-dBZ 
reflectivity  contour,  coordinates  of  the  weighted  centroids  (center  of  mass),  maxi¬ 
mum  reflectivity,  and  t lie  time  of  observation.  For  each  cell,  data  are  provided 
defining  tne  maximum  reflectivity,  ueigat  of  tne  cell's  tallest  slice,  and  the  cell' 
speed  and  direction. 

These  acquired  data  are  used  by  tne  nail  algorithm  to  determine  an  assignme 
of  one  of  three  labels  for  each  nail  indicator.  A  "V"  is  assigned  if  the  conditions 
of  tne  indicator  are  met;  an  "N"  is  assigned  if  tne  conditions  are  not  met;  finally, 
a  "L”  is  assigned  if  the  conditions  of  tae  indicator  cannot  be  tested  because  of 
insufficient  data.  An  insufficient  data  situation  occurs  when  the  predetermined 
volumetric  scan  sequence  does  not  permit  measurements  above  a  required  altitude. 

A  more  detailed  description  of  tne  hail  algorithm  is  included  in  Appendix  A. 

4.  TF.ST  CONSIDERATIONS 

The  hail  algorithm  was  tested  on  archived  radar  data  obtained  by  a  10-cm  radar 
in  Oklahoma  during  the  Joint  Doppler  Operational  Project  (JDOP).  ^  Our  test  re¬ 
sults  were  verified  bv  ground  trutli  data  provided  by  S.  Nelson  of  the  National  Severe 
Storm  Laboratory  (NSSL)  from  a  ground  observer  network  located  west  of  the  radar 
site.  We  tested  data  from  nine  selected  days  of  diverse  weather  situation  including 
a  total  of  206  observations.  A  valid  observation  required  at  least  one  of  the  two 
following  conditions:  A  radar  cell  was  identified  within  coverage  of  the  ground 
observer  network  and/or  a  ground  observer  report  was  submitted  of  precipitation 
occurring  during  tne  period  of  radar  operation. 

For  each  observation,  we  performed  our  nine  hail  indicator  tests.  The  out¬ 
come  of  each  test  and  the  associated  ground  report  for  the  observation  was  used  to 
classify  the  test  into  one  of  five  categories:  x,  y,  z,  w,  and  u.  Category  x 
(successes),  represented  tests  where  the  conditions  of  the  indicator  were  met  and 
hail  was  verified  within  an  area  defined  by  tne  30-dBZ  contour  of  the  cell.  The  time 
constraint  for  hail  verification  required  that  the  time  of  the  reported  event  be  within 
the  period  required  to  complete  a  volumetric  scan  which  was  typically  5  min.  Time 
uncertainties  reported  bv  the  ground  observers  were  added  to  the  time  constraints. 
These  uncertainties  were  generally  less  than  5  min.  Category  y  (failures )  represented 

5.  JDOP  Staff  (1979)  Final  Report  on  the  Joint  Doppler  Operational  Project  (JDOP) 

1976- 1978,  NOAA  Tech.  Memo,  EHl  NSSL-66  Norman,  OK,  84  pp. 


tests  where  the  conditions  of  tne  indicator  were  not  met  altnough  hail  was  verified. 
Category  z  (false  alarms),  represented  tests  where  the  conditions  of  the  indicator 
were  met  but  hail  was  rot  verified.  Category  w  represented  tests  where  the  condi¬ 
tions  of  the  test  were  not  met  and  hail  not  verified.  Category  u  (unknowns), 
represented  tests  that  could  not  be  made  owing  to  insufficient  data. 

We  used  the  results  of  the  tests  to  establisn  weighting  functions  for  each  of  tne 
nine  hail  indicators.  Upon  finding  these  functions,  we  were  able  to  total  a  score 
f or  each  observation  that  was  a  sum  of  the  successful  weignted  indicators.  We  then 
used  these  scores  to  identify  a  given  radar  observation  in  one  of  four  ways:  a 
hailer,  a  probable  hailer,  a  nonhailer,  or  an  unknown. 


5.  WEIGHT  DETERMINATION 

We  based  our  assigned  weighting  functions  on  ttie  critical  success  index  (CSl)^ 
for  each  indicator.  Donaldson  defined  the  CSi  as  the  ratio  of  successful  predictions 
of  a  critical  event  to  the  sum  of  successful  predictions  plus  unsuccessful  predic¬ 
tions  of  both  types.  In  terms  of  our  test  categories,  we  can  express  this  as: 

CSI  =  x/(x+ y+z)  .  (1) 

The  results  of  our  tests  were  used  also  to  determine  probabilities  of  detection 
(PODs)  and  false  alarm  ratios  (FARs). 

The  probability  of  detection  is  the  proportion  of  nail  events  correctly  predicted 
by  the  indicator  test: 

POD  =  x/ (x+y )  .  (2) 

The  false  alarm  ratio  is  the  proportion  of  false  indicator  Predictions  of  hail: 

FAR  =  z/(x+z)  ,  (3) 

A  summary  of  the  indicator  tests  and  the  calculated  success  indices  are  shown  in 
Table  1. 

It  can  be  noted  that  test  categories  w  and  u  are  not  used  to  determine  the  success 
indices.  While  the  rejection  of  u  is  obvious,  the  rejection  of  w  requires  some 
clarification.  It  would  appear  that  a  prediction  of  no  hail  when  in  fact  hail  was  not 


6.  Donaldson,  R.  J.  ,  IJver,  R.  M.  ,  and  Kraus,  IU.  (1975)  An  Objective  Evaluator 
of  Techniques  for  Predicting  Severe  Weather  Events,  Preprints,  9th  Confer¬ 
ence  on  SevereKqcal  SRmni",  Boston,  Am.  Meteorol.  Sue.,  321  -325. 
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reported,  constituted  a  success.  The  problem  lies  in  the  uncertainties  in  observing 
and  reporting  a  nail  event.  Light  hail  can  be  easily  masked  when  embedded  in  a 
heavy  shower,  or  its  occurrence  can  be  missed  if  the  edge  of  the  nail  falls  a  short 
distance  away  from  an  observer.  In  Table  1,  we  can  see  that  w  is  large  when  com¬ 
pared  with  x  +  y  +  i.  Uncertainties  in  w  would  overwhelm  the  meaning  of  any 
measure  of  the  reliability  of  the  tiail  iden*ification  technique  which  depends  on  w. 

Another  reason  to  ignore  w  is  that  the  number  of  nonhail  events  greatly  exceeds 
the  nun  ber  of  hail  events.  One  can  envision  the  use  of  a  very  insensitive  hail 
indicator  to  predict  hail.  If  w  were  counted  as  a  success,  the  overwhelming 
number  ol  w's  relative  to  a  few  failures  to  predict  hail  would  result  in  a  deceptive 
high  LSI.  Finally,  any  index  involving  w,  even  if  w  is  measured  accurately,  would 
tell  us  more  about  the  climatological  occurrence  of  hail  rather  than  the  success  of 
the  hail  identification  method. 


Table  1. 
FA  Its  for 

Summit  rv 
206  ( >bst* 

of  tiie  Indicator 
•vat  ions 

Tests 

and  tiie 

Calculated  CSIs 

PODs 

and 

Indicator 

i 

2 

.3 

4 

5 

6 

7 

8 

9 

X 

it; 

30 

10 

47 

16 

26 

8 

34 

20 

V 

l 

2 

1 

21 

1 1 

0 

5 

19 

/. 

5 

8 

1 

3.' 

8 

10 

2 

7 

25 

\v 

101 

92 

88 

113 

99 

70 

79 

87 

25 

li 

(>  ,1 

74 

106 

0 

62 

89 

1  17 

73 

117 

LSI 

0.  85  I 

0.  750 

0.  833 

0.  505 

0.  356 

0.  553 

0.  800 

0.  139 

0.  312 

poo 

0.  073 

0.  938 

0.  909 

0.  839 

0.  432 

0.  703 

1.  00 

0.  872 

0.  513 

F  Alt 

0.  122 

0.  210 

0.  091 

0.  440 

0.  333 

0.  27  8 

0.  200 

0.  17  1 

0.  556 

Although  we  recognize  that  the  hail  indicators  are  not  totally  independent 
parameters,  they  were  treated  as  such  by  normalizing  the  CSIs  to  compute  a  weight¬ 
ing  function  \\  ^  for  each  indicator  i  as  shown  by  the  expression: 


csr 


£ 


<  'SI 


(4) 


where  N  =  the  number  of  indicators. 

I  nis  equation  gives  tile  highest  weights  to  those  indicators  with  the  highest 
(  Ms  where  the  sum  of  all  the  indicator  weights  U  j  to  \V  =  1.  It  should  be  noted 
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that  this  technique  dqes  not  maximize  the  PODs  ■  -  minimize  the  FAHs.  The  t'Sl 
is  a  parameter  measuring  the  over-all  effective  ess  of  our  prediction  and  includes 
both  the  hail  events  that  we  failed  to  predict  and  our  false  alarms.  It  lias  a  range 
of  values  from  0  to  1  with  1  being  a  perfect  predictor  and  0  meaning  we  haven't 
improved  our  prediction  at  all. 

Listed  in  Table  2  are  two  sets  of  weighting  functions  for  each  hail  indicator. 
The  first  set  of  functions  was  calculated  from  the  C'Sls  determined  in  our  tests. 
From  these  results  we  can  see  that  indicators  5  and  9  are  not  useful  as  hail  pre¬ 
dictors  because  of  their  low  t'Sls.  A  second  set  of  functions  and  the  ones  we  finally 
used  were  calculated  without  indicators  5  and  9. 


Table  2.  Weighting  Functions  Determined  From  t'Sls  for  Each  Indicator 


Indicator 

12  5  4  5  6 

8  9 

Set  1  0.19  0.14  0.  18  0.06  0.  0H  0.  OH  0.16  0.  14  0.02 

Set  2  0.20  0.15  0.  18  0.  07  -  0.  OH  0.17  0.15  - 

6.  TEST  EVALUATION 

These  initial  results  show  the  relative  validity  of  the  nine  hail  indicators.  An 
important  comparison  is  one  that  we  can  make  between  the  combined  c'Sl  of  the  first 
three  indicators  anil  the  combined  C'Sl  of  the  parent  WSK-57  hail  criteria.  We  ob¬ 
tained  a  C'Sl  of  0.  8  El  compared  to  a  C'Sl  of  0.700  obtained  bv  Lemon.  This  appar¬ 
ent  discrepancy  can  be  attributed  to  the  fact  that  two  different  criteria  were  used 
to  verify  the  successful  hail  predictions.  In  the  WKS-57  hail  criteria  evaluation 
tests.  *  a  successful  prediction  required  hail  falls  with  stone  diameters  '  1.9  cm. 

No  minimum  size  limitations  were  imposed  to  verify  hail  falls  that  were  predicted 
by  our  hail  algorithm  which  would  in  effect  result  in  fewer  false  alarms  and  a  better 
C'Sl  value. 

One  of  the  unexpected  results  of  the  tests  was  the  low  effectiveness  of  the 

55-dB/.  reflectivity  (indicator  4)  as  a  dependable  hail  indicator.  In  another  study 

4 

by  Dennis,  in  Nebraska,  high  reflectivity  provided  an  excellent  indication  of  hail. 
We  can  attribute  these  conflicting  results  at  least  in  part  to  climatological  factors. 

g 

as  suggested  bv  Donaldson.  Our  data  included  many  thunderstorms  that  produced 
heavy  showers  but  no  hail.  High  reflectivities  detected  in  these  storms  resulted  in 
false  alarms  that  lowered  the  magnitude  of  the  C'Sl.  It  should  be  noted,  however, 
that  the  magnitude  of  the  I’OD  is  still  at  an  acceptable  level  and  we  should  not 
underestimate  the  usefulness  of  this  indicator. 


Another  result  in  question  is  the  unrealistic  1*01)  value  of  1.00  shown  under 
indicator  7  CIO  dH/.,  8  km).  Unfortunatelv  we  did  not  obtain  a  large  enough  data 
-ample  to  permit  a  completely  reliable  evaluation  for  this  indicator.  A  shift  of 
tust  one  sample  from  a  success  to  failure  category  would  have  had  a  large  influence 
a.  the  I'OI )  value.  Our  data  sample  was  small  because  the  heights  of  more  than  half 
>:  the  total  storm-  could  not  be  determined  owing  to  antenna  elevation  limits  during 
iota  acquisition.  \\  e  still  feel,  however,  that  this  indicator  is  an  excellent  hail 
ore  I i ■  tor. 

Indicator  6  (the  direction  of  tilt  lies  to  the  right  of  and/or  behind  the  direction 
Hnvard  which  the  cell  is  moving)  and  indicator  8  (the  direction  of  tilt  is  towards  the 
south)  are  related  in  that  tliev  both  test  for  the  direction  of  storm  tilt.  Surprisinglv, 
the  tests  show  that  indicator  8  is  significant Iv  better  than  indicator  6.  In  retro¬ 
spect  it  appears  that  the  conditions  for  testing  indicator  6  were  too  restrictive, 
narti'  Ularlv  in  determining  a  tilt  to  the  rear  of  the  storm  direction. 

Of  the  nine  indicators  tested,  5  and  !•  have  l 'Sis  below  0.500.  We  could  not. 

>f  course,  use  these  indicators  m  making  the  final  hail  determinations.  However, 
thev  still  remain  in  the  hail  algorithm  for  future  test  purposes. 

II All.  IDKM'IKICAIION 

V\  e  used  the  results  of  tin1  indicator  tests  lo  formulate  threshold  values  to 
identify  a  given  radar  observation.  Our  thresholds  were  determined  bv  a  scoring 
method  described  below. 


fable  :l.  Truth  Table  Showing  the  Results  of  the  Indicator  Tests  for  Four  Sample 
Rada r  ( >bse rvat  i o ns 


Indicator 

1 

2 

8 

4 

ft 

6 

7 

YSUM 

NST  M 

OF 

SCORK 

obs  m 

eight 

TO 

la 

18 

7 

8 

17 

15 

68 

84 

Ohs  1 

Y 

Y 

Y 

N 

T 

N 

Y 

68 

24 

92 

74 

HAIL 

(  )bs  2 

Y 

N 

.\ 

Y 

N 

N 

l 

27 

58 

85 

82 

Ohs  :< 

Y 

l 

Y 

N 

l 

1 

18 

8 

50 

84 

I’ROB 

Otis  4 

( 

r 

V 

N 

Y 

1 

T 

8 

7 

15 

58 

I'NKN 

In  Table  8.  VSTM  is  the  weighted  sum  of  tests  that  have  met  their  criteria  (Y)s; 
NSt  IW  is  the  weighted  sum  of  tests  that  have  not  met  their  criteria  (NIs.  The  (FIs 


are  tests  that  could  not  be  performed  because  of  insufficient  data;  CF  is  the  sum  of 
YSl'M  and  NSU1V1.  Since  CF  is  the  maximum  possible  weighted  sum  for  tests  that 
could  be  performed  for  a  given  cell,  we  can  consider  CF  as  a  Confidence  Factor 
for  that  cell;  CF  =  100  if  all  of  the  tests  are  performed.  A  relative  score  (SCORE) 
is  determined  by  a  ratio  of  the  sum  of  the  successful  tests  to  the  total  possible  sum 
or 


SCORE  -  (YSUM/CF)  X  100  .  (5) 

After  finding  the  SCORES  for  each  of  the  206  observations  in  our  tests,  we 
endeavored  to  find  the  best  SCORE  thresholds  that  could  identify  storms  with  the 
highest  probability  of  success.  The  CSIs,  PODs,  and  FARs  were  calculated  for 
ten  ranges  of  the  SCORE  and  for  three  ranges  of  the  confidence  factor  (CF)  as  shown 
in  Table  4.  Once  again,  for  each  SCORE  and  CF  range  we  define  x  as  cells  in  which 
hail  indicator  conditions  were  met  and  hail  confirmed;  y  as  cells  in  which  hail  indi¬ 
cator  conditions  were  not  met  but  hail  reported;  z  as  cells  in  which  the  hail  condi¬ 
tions  were  met  but  hail  not  confirmed;  u  as  cells  in  which  the  hail  indicator  condi¬ 
tions  could  not  be  tested;  and  w  as  cells  in  which  the  indicator  conditions  were  not 
met  and  hail  not  reported.  The  confidence  level  ranges  were  arbitrarily  chosen  to 
show  the  trend  of  success  indices  as  the  confidence  level  increases. 

A  comparison  of  the  CSIs  vs  SCORE  for  the  three  CFs  is  shown  by  the  graph  in 
Figure  3.  From  the  graph,  we  can  see  that  the  best  threshold  for  identifying  hail¬ 
storms  for  all  three  CFs  is  a  SCORE  that  is  at  least  equal  to  60.  We  can  also  see 
that  as  the  CF  becomes  greater  we  increase  our  CSI  but  lose  some  of  our  cell  popu¬ 
lation.  A  good  compromise  for  identifying  hailstorms  turns  out  to  be  a  SCORE  2  60 
with  a  CF  >  50. 

Figure  4  provides  a  flow  diagram  showing  the  required  conditions  for  deter¬ 
mining  our  four  storm  identifiers.  Radar  observations  in  which  only  25  percent  or 
less  of  weighted  tests  could  be  made  (CF  -  25)  are  disregarded  and  labeled  as 
unknown  (UNKN).  Remaining  storms  with  scores  -  60  are  labeled  as  hailer  (HAIL) 
if  CF  >  50  or  as  probable  hailers  (PROB)  if  CF  S  50.  Observations  with  scores  <  60 
are  considered  to  be  nonhailers.  Going  back  and  applying  the  hail  algorithm  to  our 
data  that  included  206  radar  observations,  we  obtained  the  results  shown  in 
Table  5. 
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Figure  4.  Conditional  Flow  Diagram  for  Identifying 
Hailstorms 
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Table  4.  Hail  Success  Indices  vs  SCORES  for  Three  Confidence  Factor  (CF) 
Ranges 


CF  >  0 

SCORE 

;  o 

2  10 

2  20 

2  30 

2  40 

2  50 

2  60 

2  70 

2  80 

2  90 

X 

56 

53 

53 

52 

51 

51 

51 

48 

42 

39 

y 

0 

3 

3 

4 

5 

5 

5 

8 

14 

17 
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150 

42 

37 

35 

32 

28 

24 

23 

23 

21 

u 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

w 

0 

108 

1  13 

115 

118 

122 

125 

127 

127 

129 

CSI 

0.  27 

0.  54 

0.  57 

0.  57 

0.  58 

0.  61 

0.  64 

0.  61 

0.  53 

0.  51 

POD 

1.  00 

0.  95 

0.  95 

0.  93 

0.  91 

0.  91 

0.  91 

0.  86 

0.  75 

0.70 

FAR 

0.  73 

0.  44 

0.  41 

0.40 

0.  39 

0.  35 

0.  32 

0.  32 

0.  35 

0.  35 

CF  >  25 

SCORE 

5  0 

B9 

>  20 

2  30 

2  40 

2  50 

2  60 

2  70 

2  80 

2  90 

X 

37 

37 

37 

37 

35 

35 

35 

32 

26 

23 

y 

0 
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0 

2 

2 

2 

5 

1  1 

14 

Z 

106 

21 

17 

15 

12 

8 

5 

3 

2 

2 

u 

63 

63 

63 

63 

63 

63 

63 

63 

63 

63 

w 

0 

85 

89 

9  1 

94 

98 

101 

103 

104 

104 

CSI 

0.  26 

0.  64 

0.  69 

0.7  1 

0.7  1 

0.  78 

0.  83 

0.  80 

0.  67 

0.  59 

POD 

1.  00 

1.00 

1.  00 

1.  00 

0.  95 

0.  95 

0.  95 

0.  87 

0.  7  0 

0.  62 

FAR 

0.74 

0.  36 

0.  32 

0.  29 

0.  26 

0.  19 

0.  12 

0.  09 

0.  07 

0.  08 

C  F  >  50 

SCORE 

>  0 

?  10 

>  20 

2  30 

2  40 

2  50 

?  60 

-70 

■  80 

'  90 

X 

33 

33 

33 

33 

31 

31 

3  1 

28 

23 

20 

y 
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0 

2 

2 

2 

5 

10 

13 

Z 

97 

13 

10 

8 

5 

3 

2 
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0 

0 

u 

76 

76 

76 

76 

76 

76 

76 

76 

76 

76 

w 

0 

84 

87 

89 

92 

94 

95 

97 

97 

97 

CSI 

0.  25 

Ksa 

0.  86 

0.  89 

0.  85 

0.70 

0.  61 

POD 

1. 00 

0.  94 

0.  94 

0.  94 

0.  85 

0.  70 

0.  61 

FAR 

0.  75 

m 

0.  29 

0.  14 

0.  09 

0.  06 

0.  00 

0.  00 

0.  00 
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Table  5.  Hail-  Identification  Summary  and  Success  Indices  for  206  Radar 
Observations 


Cell 

Identifiers 

Number  of 
Observations 
Identified 

Hail 

Occurrences 

Failed  to 
Detect  Hail 

False  Alarms 
of  Hail 

Hailers 

33 

31 

2 

2 

Prob  Hailers 

7 

4 

0 

3 

Nonhailers 

103 

2 

Unknowns 

63 

19 

Cell 

Identifiers 

Successes 

(x) 

Failures 

(y> 

False 

Alarms 

(z)  POD 

FAR 

CSI 

Hailers 

31 

2 

2 

0.  939 

0.  061 

0.  886 

Hailers  and 
Prob  Hailers 

35 

2 

5 

0.  946 

0.  125 

0.  833 

Out  of  206  observations,  the  hail  algorithm  identified  33  cells  as  HA1LERS. 

Only  two  of  the  identified  HAILERS  were  false  alarms  and  only  two  hail  occurrences 
were  not  identified  by  the  algorithm  to  give  as  a  probability  of  detection  (POD  =  0.  939, 
a  false  alarm  ratio  (FAR)  =  0.061  and  a  critical  success  index  (CS1)  =  0.886.  We 
were  able  to  identify  an  additional  seven  PROBABLE  HAILERS  by  using  a  lower 
confidence  factor  but  we  also  increased  our  false  alarms.  When  using  the  lower 
confidence  factor  we  in  effect  improved  our  probability  of  detection  to  0.  946  but 
increased  our  false  alarm  ratio  to  0.  125  and  decreased  our  critical  success  index 
slightly  to  0.  833.  Out  of  the  103  observations  identified  as  NONHAILERS,  only  two 
were  reported  as  hailers.  We  did  not  use  these  observations  to  calculate  our 
success  indices  because  of  the  uncertainties  in  the  reporting  of  nonsevere  events 
but  the  correct  identification  of  a  NONHAILER  is  in  Tact  a  successful  prediction. 

All  of  the  results  of  the  most  current  analysis  are  stored  on  disc  for  subsequent 
call  by  the  MRAS  print  and  plot  modules.  Our  print  output  lists  the  results  of  the 
analysis  in  a  format  similar  to  the  one  shown  in  Table  3.  The  plot  output  is  dis- 
displayed  on  a  CRT  monitor  in  conjunction  with  the  track  output.  In  addition  to 
other  track  information,  each  cell  is  annotated  with  a  symbol  H,  P.  U  or  a  blank 
to  indicate  if  the  cell  was  identified  as  a  hailer.  a  probable  bailer,  an  unknown  or 
a  nonhailer,  respectively,  and  the  associated  score. 


8.  OOiNCl.l  SIGNS 


We  have  shown  that  hailstorms  can  be  automatically  identified  with  a  high  degree 
>f  reliability  by  computer  processing  of  real-time  radar  data.  The  advantages  over 
the  past  hail  identification  techniques  are  twofold.  First,  the  computer  processing 
allows  testing  of  many  hail  criteria.  We  were  able  to  make  at  least  nine  tests 
without  timing  limitations.  Although  we  did  exclude  two  of  our  initial  tests,  the 
results  indicate  that  our  <_'SI  increases  as  a  function  of  the  number  of  valid  tests. 

Our  second  advantage  lies  in  the  nature  of  the  automatic  technique,  allowing  an 
analysis  of  every  storm  within  the  radars  coverage  area. 

The  only  apparent  limitations  of  the  automatic  hail  detection  technique  is  the 
data  acquisition  constraint  by  the  elevation  scan  sequence.  A  typical  scan  sequence 
used  to  acquire  the  data  for  our  study  was  composed  of  five  elevations  that  included 
0.  2,  0.  5,  1.  5,  3.  0.  and  5.  0  degrees.  At  a  maximum  elevation  of  5.  0  degrees,  we 
could  not  measure  tops  of  storms  above  8  km  when  the  ground  ranges  of  the  storms 
were  less  than  86  km.  We  were  able  io  test  storms  within  this  range  with  other 
hail  indicator  tests  but  with  less  confidence. 

In  this  study  we  were  able  to  correlate  storms  producing  hail  without  size  dis¬ 
tinction  to  a  SCORE  threshold  resulting  from  successful  indicator  tests.  A  natural 
extension  of  this  work  would  be  to  correlate  also  hail  size  to  the  magnitude  of  the 
SCORE.  Recalling  that  the  SCORE  is  actually  a  ratio  of  the  positive  tests  to  the  total 
tests  performed,  we  can  use  only  SCORES  where  100  percent  of  the  tests  were  per- 
formedto  make  a  valid  SCORE  vs  hail  size  correlation.  In  our  study  we  found  93  obser¬ 
vations  in  which  100  percent  of  thetests  were  performed  but  only  6  of  these  observa¬ 
tions  were  hailers.  We  tried  to  correlate  the  SCORE  of  these  observations  to  hail 
size  but  the  results  were  inconclusive.  Our  inability  to  make  the  SCORE  vs  hail 
size  correlation  can  be  linked  to  our  inability  to  make  all  of  the  indicator  tests 
particularly  in  cells  within  86  km. 

It  should  be  pointed  out  that  the  elevation  scan  sequence  used  to  acquire  our 
data  has  been  designed  with  other  objectives  in  mind.  We  intend  to  continue  our 
hail  studies  using  a  tailored  elevation  scan  sequence  that  will  enable  complete  testing 
of  more  cells.  One  of  the  primary  objectives  of  future  studies  will  be  to  extend 
>ur  hail  identification  capability  to  include  a  prediction  of  hail  size. 
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Appendix  A 

Hail  Algorithm  Description 


Al.  FUNCTIONAL  DESCRIPTION 

The  purpose  of  this  algorithm  is  to  predict  one  of  the  following  four  cases: 

•  a  given  storm  produces  hail, 

•  a  given  storm  probably  produces  hail 

•  a  given  storm  does  not  produce  hail, 

•  a  given  storm  cannot  be  analyzed  due  to  lack  of  insufficient  data. 

This  algorithm  analyzes  storm  data  available  in  a  specific  format  as  described 
by  the  following.  A  storm  is  defined  as  a  three-dimensional  region  of  significant 
reflectivity  values  (2  30  dBZ).  It  is  assumed  to  be  made  up  of  two  dimensional 
circular  storm  components  occurring  at  different  elevation  angles  of  radar  obser¬ 
vation.  A  storm  component  has  associated  with  it,  a  centroid,  a  maximum  re¬ 
flectivity  value  and  an  area.  The  centroid  is  represented  in  a  Cartesian  coordinate 
system  with  the  radar  at  the  origin.  The  X-axis  denotes  east-west  directions  and 
the  Y-axis  denotes  north-south  directions.  The  maximum  reflectivity  value  in  a 
storm  component  is  defined  as  the  largest  value  of  reflectivity  obtained  from  the 
radar  for  the  resolution  volumes  within  the  identified  storm  component.  The  area 
of  a  storm  component  is  defined  as  the  area  of  the  circle  representing  the  actual 
storm  component.  The  centroid,  the  maximum  reflectivity  and  the  area  for  each 
storm  component  within  a  storm  are  obtained  (refer  to  STORM  STRUCTURE, 
NX-DR-03-009)  as  inputs  to  this  algorithm.  Additional  storm  parameters  will  be 


23 


used  as  prime  inputs.  These  are  the  estimated  speed  and  direction  of  storm  move¬ 
ment,  the  maximum  reflectivity  value  within  a  whole  storm,  and  the  height  of  the 
storm's  highest  detectable  storm  component. 

Two  definitions  are  necessary  to  understand  the  process  of  analyzing  the 
storm  data.  First  is  the  overhang.  A  storm  is  said  to  have  an  overhang  if  the  edge 
of  the  storm  component  at  height  between  5  km  and  12  km  extends  beyond  the  edge 
of  a  storm  component  at  the  lowest  elevation  by  at  least  4  km.  Second,  a  storm  is 
said  to  have  a  tilt,  if  the  centroid  of  a  midlevel  storm  component  is  to  the  right  or 
to  the  rear  of  the  lowest  level  storm  component. 

Eight  hail  predictors  are  used  to  determine  whether  a  given  storm  produces 
hail,  will  probably  produce  hail,  or  will  not  produce  hail.  A  weight  is  associated 
with  each  predictor  based  on  the  empirical  study  of  storms  at  AFGL.  The  pre¬ 
dictors  and  their  weights  are: 

1.  The  highest  detectable  storm  component  is  at  least  8  km.  (Weight  =  171 

2.  The  maximum  reflectivity  value  within  a  storm  is  greater  than  55  dBZ. 
(Weight  =  7) 

3.  The  centroid  of  the  lowest  level  storm  component  is  to  the  north  of 
a  centroid  at  any  higher  level  (Weight  =  15) 

4.  The  direction  of  tilt  is  between  45  and  100  degrees  to  the  right  of 
the  direction  of  storm  movement.  (Weight  =  8) 

5.  The  area  of  a  storm  component  at  any  midlevel  elevation  is 
greater  than  that  at  the  lowest  elevation.  (Weight  -  0) 

6.  The  maximum  reflectivity  in  a  storm  component  at  midlevels 
(that  is,  5  to  12  km  heights)  is  at  least  50  dBZ.  (Weight  =  20) 

7.  The  overhang  at  the  midlevel  extends  to  at  least  4  km  beyond  a 
lower  level  storm  component.  (Weight  =  15) 

8.  The  highest  elevation  storm  component  exists  above  a  midlevel 
overhang.  (Weight  =  18) 

The  above  predictors  can  be  identified  in  two  ways.  A  positive  identification 
occurs  when  sufficient  data  exists  to  test  for  the  presence  of  the  corresponding  con¬ 
dition  associated  with  the  predictor.  Second,  a  probable  identification  occurs  when 
sufficient  data  does  not  exist  (for  example,  the  radar  could  not  take  measurements 
beyond  a  fixed  elevation  angle).  The  algorithm  tests  for  tnese  two  types  of  pre¬ 
dictors  and  then  computes  two  sums,  one  using  only  positive  predictors  and  another 
using  only  probable  predictors.  The  decision  of  labelling  a  storm  as  positively, 
probably,  or  not  hail  producing,  or  storm  data  is  insufficient  for  hail  analysis  is 
then  made  as  indicated  by  the  procedure  section. 
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Al.l  Source 


This  algorithm  has  been  implemented  by  Air  Force  personnel  at  the  Air  Force 
Geophysics  Laboratory  (AFGL),  Sudbury,  Massachusetts 

A1.2  Procesaing  Environment 

The  HAIL  algorithm  is  implemented  as  one  of  four  special  analysis  modules 
that  operate  under  the  Modular  Radar  Analysis  Software  Systems  (MRASS)  (see 
Reference  2  in  Section  1)  implemented  at  AFGL,  Sudbury.  The  algorithm  requires 
outputs  from  the  modules  implementing  algorithms  described  in  STORM  STRUCTURE 
(NX -DR- 03 -009),  STORM  FORECAST  (NX-DR-03-008),  and  STORM  CENTROIDS 
(NX -DR- 03 -005). 


A2.  INPUTS 


A2.1  Identification 
AREA 

DIRECTION 


ELEVATIONS 

MAXIMUM  STORM 
REFLECTIVITY 
SPEED 


STORMS 


STORM  COMPONENT 
REFLECTIVITY 
STORM  TOP 


I 

; 


Areas  of  the  storm  components  of  an  identified 
as  given  by  STORM  STRUCTURE  (NX-DR-03-009). 
Direction  of  motion  of  each  identified  storm  as 
given  by  STORM  FORECAST  (NX-DR-03-008).  If 
the  direction  is  not  determined  by  STORM  FORECAST, 
then  its  value  is  set  as  unknown. 

Elevation  angles  of  each  radar  scan  constituting 
one  volume  scan. 

Maximum  detected  reflectivity  value  among 
all  storm  components  of  each  storm. 

Speed  with  which  each  identified  storm  is 
moving  as  given  by  STORM  FORECAST 
(NX-DR-03-008).  If  the  speed  is  not  determined 
by  STORM  FORECAST,  then  its  value  is  set 
as  unknown. 

Identifiers  for  three-dimensional  regions 
characterized  by  a  number  of  storm  components 
taken  at  successive  elevations  and  with 
reflectivity  values  above  a  given  threshold. 

Maximum  reflectivity  value  detected  in  an 
individual  storm  component. 

Altitude  of  the  centroid  of  the  highest  detectable 
storm  component  in  each  storm.  If  the  storm 
component  is  detected  at  the  highest  elevation 
angle  the  value  is  signed  as  a  negative  number. 
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times 


\- POSITION 


\  -I'OSIi'KlN 


W  EIGHT 


VJ .  2 

ARE A  Us  acquired  as  an  output  from  STORM  STRUCTURE  (NX-DR-03-009). 

ELE\  ATIONS  and  TIMES  are  acquired  directly  as  measured  values  from  the 
Doppler  radar. 

SPEED  and  DIRECTION  are  acquired  as  outputs  from  STORM  FORECAST 
(NX-DH-03-008). 

STORMS  are  acquired  as  outputs  from  STORM  CENTROIDS  (NX-DR-03-005). 

N- POSITION,  \  -POSITION,  STORM  TOP,  MAXIMUM  STORM  REFLECTIVITY, 
and  STORM  COMPONENT  REFLECTIVITY  are  acquired  as  outputs  from  STORM 
STRUCTURE  (NX-DR-03-009). 

WEIGHT  is  a  system  supplied  parameter  whose  value  is  based  on  the  study  of  hail 
producing  storms. 


=  Observation  times  at  the  beginning  of  each  radar 
scan,  one  for  earn  elevation  angle. 

=  X-positions  of  the  mass  weighted  centroids 
(centers  of  mass)  of  the  storm  components 
observed  at  TIMES. 

-  \  -positions  of  the  mass  weighted  centroids  (centers 
of  mass)  of  the  storm  components  observed 

at  TIMES. 

-  Weight  factor  associated  with  each  of  the 
hail  predictors  identified  in  Section  Al. 


AH.  PROCEDURE 
AH.  1  Algorithm 

BEGIN  ALGORITHM  (HAIL) 

1.0  COMPUTE  (AVERAGE  SPEED) 

2.0  COMPUTE  (AVERAGE  DIRECTION) 

3.  0  DO  FOR  ALL  (STORMS) 

3.  1  IF  (SPEED  is  unknown) 

THEN  (Set  SPEED  to  AVERAGE  SPEED) 
END  IF 
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3.  2  EE  (DIRECTION  is  unknown  UK  DIRECTION  differs  from 
tile  AVERAGE  DIRECTION  be  90°  or  more) 

THEN  (Set  DIRECTION  to  AVERAGE  DIRECTION) 

END  IE 

3.3  IE  (STORM  TOP  is  not  known) 

THEN  (Identify  probable  HAIL  PREDICTOR  #1) 

END  IE 

3.4  IP  (STORM  TOP  is  at  least  8  km) 

THEN  (Identify  positive  HAIL  PREDICTOR  #1) 

END  IE 

3.  5  IP  (MAXIMUM  STORM  REFLECTIVITY  is  greater  than 
55  dHZ) 

THEN  (Identify  positive  HAIL  PREDICTOR  #2) 

END  IF 

3.  6  COMPUTE  (N-SPEED_ 

3.7  COMPUTE  (V -SPEED) 

3.  8  DO  FOR  ALL  (ELEVATIONS) 

3.  8.  1  COMPU  TE  (RADIUS) 

3.8.2  COMPUTE  (RANGE) 

3.8.3  COMPUTE  (HEIGHT  ) 

END  DO 

3.  9  DC )  (ELEV  ATIONS)  FROM  (Second  lowest)  TO  (Highest) 

3.  9.  1  COMPUTE  (X- DISPLACEMENT) 

3.  9.  2  COMPUTE  U  -DISPLACEMENT) 

3.  9.  3  COMl'U  TL  (DELTA  ANGLE) 

3.9.4  COMPUTE  IDISIANCE) 

3.  9.5  IF  (3  -  DISPLACEMENT  is  negative) 

THEN  ‘Identify  positive  HAIL  PREDICTOR  #3) 

END  11 

3.  9.8  IF  (DELTA  ANGLE  is  between  45  and  180  degrees) 
THEN  (Identify  positive  HAIL  PREDICTOR  #4) 

END  IF 

3.9.7  IF  (RADIUS  of  current  storm  component  is  greater 
than  RADIUS  of  lowest  elevation  storm  component) 
THEN  (Identify  positive  HAIL  PREDICTOR  #5) 

END  IF 

3.  9.  8  IF  (HEIGHT  of  the  current  storm  component  is  not 

between  5  and  12  km  AND  STORM  TOP  is  not  known) 
THEN  (Identify  probable  HAIL  PREDICTOR  #6) 
(Identify  probable  HAIL  PREDICTOR  #7' 

END  IF 
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■.  :•  11  .  ill.li.ir:  "l  tin-  ,  .i!  it-!!*  storm  i  ompmient  is 

Iji  'mitii  5  i!i<!  1-  kni  ANU  ST'  >KM  COMPONUNl 
Kill  I.E<  ilNIl'i  in  erf  itei'  tli.m  50  dll/.) 

JMUN  l  Identify  positive  1 1 A 1 1 .  PREDICTOR  »  6  I 
I.MJ  II 

3.  ■*.  10  II  (HER, ill  of  tin  loucs*  elevation  storn.  component 

is  less  than  5  km  A X I )  HEIGHT  of  the  current  storm 
component  is  between  5  -ind  12  km l 
THUN  A  »i\l  I M  ill  IO\  UKHANGl 

1 1  (OVERHANG  is  ere. iter  t;i;m  -4  km> 

X 1 1 1 .  \  l  identif  \  positive  I!  All.  I'KKUIt'TuH  r7l 
I  All  II 

END  I  )■  > 

3.  10  II  ill). Kill!  of  trie  las*  storm  component  is  ibo.e  the  HEIGHT 
■  tr.i  storn.  component  ‘iia*  Ins  '.he  OVERHANG' 

1  1 1 1 :  \  CQiM'l  1  I.  'siOKM  ini'  D1>1  ANCE* 

1 1  (STORM  1  ( )1 '  DISi  AM  I.  is  less  than  or  equal  to 
‘he  RADII  a  ol  the  storm  component  that  contains 
tin-  <  >\  UK  I  LAM  U 

THUN  (identify  positive  ILAII.  PREDICTOR  ?8) 

END  11 

l  ,\  i.)  a 

H .  11  11_  (SKiKiM  ini'  is  not  known) 

1  I1UN  (Identify  probable  HAIL  PREDICTOR  »B) 

UNU  II 

3.  12  II  (STORM  TOP  is  unknown  A N U  HEIGHT  of  the  highest 
elevation  storm  component  is  less  than  5  km) 

THUN  (Identify  probable  ILAIL  PKUUICTOKS  #2,  #3,  #4,  and  f?5) 
UNU  IP 

3.  13  COMPL  IE  ( POSITIN' U  WEIGHT) 

3.  14  COiUPl  1  U  ( PROHA  HI.  U  tt  EIGHT) 

3.  15  COMP l  I  K  (CON F ID U N CU  FACTOR) 

3.  16  IF  (CONI  IDENCE  FACTOR  is  greater  than  25) 

THUN  COMPUTE  (SCORE) 

ELSE  (LABEL  the  current  STORM  as  having  insufficient  data) 
END  II 
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3.  17  IF  (SCORE  is  less  than  60) 

THEN  (LABEL  the  current  STORM  as  a  non-hail  producer) 

IF  (CONFIDENCE  FACTOR  is  greater  than  50) 

THEN  (LABEL  the  current  STORM  as  a  hail  producer) 

ELSE  (LABEL  the  current  STORM  as  a  probable  producer  of  hail) 
END  IF 
END  IF 

3.  18  VV RITE  (LABEL) 

END  DO 

END  ALGORITHM  (HAIL) 

A3.2  Computation 

A3.  2.  1  NOTATION 

3  -  AVERAGE  SI’EED  OF  N  storms  in  km/sec. 

N  =  Number  of  storms  whose  speeds  are  known. 

S.  -  Speed  of  if)l  storm  from  Sl'EED  in  km/sec. 

"D  =  AVERAGE  DIRECTION  of  N  storms  in  radians. 

D.  -  Direction  of  storm  from  DIRECTION  in  radians. 

S.V  =  X-SPEED,  X-component  of  the  speed  of  the  i111  storm  in 

km  /sec. 

S\V  =  V -SPEED,  V -component  of  the  speed  of  the  i11'  storm 

in  km /see. 

R.  =  RADIUS  of  the  storm  component  of  the  itfl  storm  at 

t.l 

elevation  in  km. 

AR.  j  =  Area  of  the  storm  component  of  the  i**1  storm  at  elevation 
in  km"  from  AREA. 

Pi  =  Mathematical  constant  having  a  value  of  3.  1416. 

RC.  .  =  RANGE  to  the  centroid  of  the  storm  component  of  the  i  1  storm 

11  ,  til  ,  4  , 

at  |  elevation  in  km. 

XC  =  X-P()SITK)N  of  the  centroid  of  the  storm  component  of 

’•*  ..  th  ,  th  ,41 

the  l  storm  at  ,|  elevation  in  km. 

YC\  j  =  Y- POSITION  of  the  centroid  of  the  storm  component  of 

the  i*^1  storm  at  elevation  in  km. 

H.  =  HEIGHT  of  Die  centroid  of  the  storm  component  of 

.  .th  th  .  .  , 

ttie  l  storm  at  i  elevation  in  km. 

( p  =  The  j  elevation  angle  in  radians  available  from 

ELEVATIONS. 


dX 

1.1 

dT. 

.1 

dV 

1.1 

dQi. 

HU  . 
i.l 

O. 

i.l 

OD. 

i 

PSW 

P  K\V 

POS 


NOTE: 


=  X-DISPLACEMEN  T.  relative  distance  in  the  X-direction 
between  the  location  of  the  centroid  of  the  i^1  storm 
at  j  elevation  and  the  location  of  the  centroid  of  the 
i**1  storm  at  lowest  elevation  in  km. 

=  Time  elapsed  between  radar  scan  at  j**1  elevation 
and  radar  scan  at  lowest  elevation  in  seconds. 

=  V -DISPLACEMENT,  relative  distance  in  the 
Y -direction  between  the  location  of  the  centroid 
of  ttie  i*n  storm  at  j  '  elevation  and  the  location  of  the 
centroid  of  the  storm  a*  lowest  elevation  in  km, 

=  DELTA  ANGLE,  the  horizontal  angle  formed  by  a 
line  from  the  centroid  of  the  storm  component  of 
the  i* '*  storm  at  /  '  elevation  and  the  centroid  of  the 
storm  component  of  the  same  storm  at  lowest  elevation 
to  a  line  that  is  parallel  to  the  direction  of  the  storm's  movement. 

-  The  horizontal  DISTANCE  between  the  centroid  of  the 
storm  component  of  the  i^  ’  storm  at  i  h  elevation  and  the 
centroid  of  the  storm  component  of  the  same  storm  at 
lowest  elevation  in  km. 

=  OVERHANG,  the  distance  between  the  edge  of  the  storm 
component  of  tin-  i^1  storm  at  i**1  elevation  and  the 
furthesl  edge  of  tile  storm  component  of  the  same  storm 
at  lowest  elevation  in  km. 

;  STOKIU  TOP  DISTANCE,  the  horizontal  distance  between 
the  centroid  of  a  storm  component  of  the  i**1  storm  at  .< 
midlevel  altitude  (5  km  to  12  km)  and  the  centroid  of  the 
storm  component  of  the  same  storm  at  the  highest  elevation, 
in  km. 

=  POSITIVE  WEIGHT,  total  weight  of  all  the  positive  hail 
predictors  identified. 

=  PROBABLE  WEIGHT,  total  weight  of  all  the  probable 
hail  predictors  identified. 

=  POSSIBILITY,  a  number  indicating  the  measure  of  the 
probability  of  labelling  a  STORM  to  be  a  potential  hail 
producer. 


This  algorithm  has  been  implemented  on  a  32-bit  minicomputer. 
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A3.  2.  2  SYMBOLIC  FORMULAS 
COMPUTE  (AVERAGE  SPEED) 


COMPUTE  (AVERAGE  DIRECTION) 


COMPUTE  (X-SPEED) 

SX^  =  (ST  (sin  D.) 

COMPUTE  (Y-SPEED) 

SY.  =  (S.)  (cos  D.) 

COMPUTE  (RADIUS) 

Rj  =  (AR../Pi)1/2 
COMPUTE  (RANGE) 

RCy-  [(SCyHXC^MYCyHYCyH1'2 
COMPUTE  (HEIGHT) 

H..  =  {(RC..)2/[<2)  (1.  21)  (637  1)]^  +  { <RC. .)  (sin  <0  .»} 

where  1.  21  is  the  index  of  refraction  and  637  1  is  the 
radius  of  earth  in  km. 


COMPUTE  (X-DISPLACEMENT) 

"ij  =  XCij  -  (XCi(lowest)  +  <SXi>  Wj» 
COMPUTE  (Y -DISPLACEMENT) 


dY..  = 


U...  +  (SY.)  (dT.)) 

i(lowest)  i  l 
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COMPUTE  (DELTA  ANGLE) 

dQjj  =  [tan"1  (cUCj/dY.j)]  -  H 

COMPUTE  (DISTANCE) 

HDjj  =  KdX..)2  +  (dY^)2)172 

COMPUTE  (OVERHANG) 

O. .  =  HD. .  +  R. .  -  R...  .. 

ij  ij  lj  l(lowest) 

COMPUTE  (STORM  TOP  DISTANCE) 

OD.  =  [(dX.,  ...  -  dX...  .  ,  ,.)2  +  (dY.,  ...  -  dY...  .  .  „.u)2]1/2 

1  1  i(mid)  i(highest)  i(mid)  i(highest))  J 

COMPUTE  (POSITIVE  WEIGHT) 

PSW  =  Sum  of  the  WEIGHT  of  positive  HAIL  PREDICTOR  #i, 
where  i  is  index  for  positive  predictors 

COMPUTE  (PROBABLE  WEIGHT) 

PRW  =  Sum  of  the  WEIGHT  of  probable  HAIL  PREDICTOR  #j. 
where  j  is  the  index  for  probable  predictors 

COMPUTE  (CONFIDENCE  FACTOR) 

CF  =  100  -  (PROBABLE  WEIGHT) 

COMPUTE  (SCORE) 

SCR  =  ((POSITIVE  WEIGHT)/(CONFIDENCE  FACTOR))*lOO 

A4.  INFERENCES 

A4.1  limitation* 

This  algorithm  is  limited  to  three-dimensional  reflectivity  structure  and 
assumes  that  ail  storm  component  are  circular.  To  provide  optimum  results,  the 
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